from scipy.integrate import quad

# Define a function to integrate
def integrand(x):
    return np.exp(-x**2)

# Compute the integral
result, error = quad(integrand, -np.inf, np.inf)
print("Integral value:", result)